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Abstract 

The partition function of the 2D Ising model coupled to an external 
magnetic field is studied. We show that the sum over the spin variables 
can be reduced to an integration over a finite number of variables. This 
integration must be performed numerically. But in order to reduce 
the partition function we must introduce as many different coupling 
constants as spin variables. The total memory that we need in order 
to store these coupling constants imposed important restrictions on 
the number of spin variables. 
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1. Introduction 

Physical models in Statistical Mechanics and Quantum Field The- 
ory are defined by integrals or sums over an infinite number of variables 
like: 

lim ^.-Jap^i,-,^}. (1) 

171 <7jV 

The weight W is the energy in Statistical Mechanics and the classi- 
cal action in Quantum Field Theory. The variables {a{} are coupled 
through terms like 



E 



0"icri+i- (2) 



In general we cannot solve exactly (^Q). Actually these sums can 
be solved exactly only if they can be transformed into sums over un- 
coupled variables. For instance, in free field models the Fourier coef- 
ficients of the field variables are not coupled. In Statistical Mechanics 
the ID Ising model^P and the 2D Ising model without magnetic field 
are equivalent to free fermion models [2J El ■ 

In this paper we are going to show that we can decouple the spin 
variables of the partition function of the 2D Ising model with an ho- 
mogeneous magnetic field. This is an interesting result because the 
2D Ising model with magnetic field cannot be solved exactly. 

In order to decouple the spin variables we are going to use the 
following trick. Let us remark that if the interacting term between 
the variables is given by 




(3) 

then we can decouple the variables using the following identity 




ocexp^ VaJ }. (4) 



Let us consider 

t 1 ™ h £ dUJ (? e ^ %i ) (? e ~ tuj6k ' lak ) ' (5) 

where {6j} is a set of different constants. In other words 

Oi^Oj if i + j- (6) 
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Let us remark that 

lim A frf^Vo (7) 
only if 6 is zero. Hence © becomes 

^JX^A-i)^. (8) 

j k 

But the constants {9i} are different. Hence the above expression be- 
comes 

X X hk-i^k-l = X W+l- ( 9 ) 
Hence we can represent the interacting term 

X a i a i+i ( 10 ) 

i 

as a double sum over independent indices and we can use @ in order 
to decouple the spin variables. 

In the next section we will apply this trick to the 2D Ising model 
with a magnetic field. 

Once the spin variables are decoupled we can integrate over all 
the spin variables. Hence we will show that the partition function of 
the 2D Ising model is given by an integral over a finite number of 
variables. 

From a theoretical point of view this result is very interesting be- 
cause the 2D Ising model with magnetic field cannot be solved exactly 
but we will show that the sum over the arbitrary number of spin vari- 
ables can be reduced to an integral over a finite number of variables. 

But from a practical point of view the trick developed in this paper 
has an important drawback. We will show that the set of coupling 
constants {9i} must take their values on a very wide range of numbers. 
Each coupling constant Oi has a fixed value but if the number of spin 
variables is large then some of the coupling constants must take very 
large values. Therefore the number of spin variables cannot be large. 

Therefore this approach can only be used for small lattices. But 
with the hypothesis of finite-size scalingjl] we do not need very large 
lattice in order to study the physics of the model at the thermo- 
dynamic limit. Moreover for some boundary conditions the model 
reaches the thermodynamic limit faster HJ. 
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In[H| it has been shown that the dimer approach [S] can be used 
in numerical calculation for small lattices. But the dimer approach 
cannot be used when a magnetic field is present or in 3D. Actually 
the approach that we study in this paper can be generalized to the 3D 
Ising model. 

2. The partition function of the 2D Ising model 

Let us consider the set {<7j} of spin variables: 

<>i = ±1 (11) 

defined over vertices of the square lattice. Index i labels columns and 
j labels rows. Let us define the energy: 

E = P<*ij(<Tij+l + + acr i,j- (I 2 ) 

The 2D Ising model with coupled to an external magnetic field is 
defined by the partition function: 

Z = 5>xp[£]. (13) 

Because the spin variables are defined by (fTTj) the partition function 
is also given by 

Z = (cosh p)*** (cosh a) Ny Y^ H 

(1 + zaijaij + i)(l + zai tj a i+ ij)(l + ha itj ), (14) 

where Nl is the number of links and Ny is the number vertices of the 
lattice. And 

z = tanh P 

h = tanh a. (15) 

For each o"jj let us define four new spin variables: afj, afj, af-j and 
a\y Now let us define the partition function: 

Z =AY,aUi,j (l + V^X^U + V^M^) 
(1 + y/zO-ijCTi^il + sfzOijofj) 

(1 + ha itj )(l + a^af +1>j )(l + rr/>^ ; . , ). (16) 
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where 

/ 1 \ 4TV V 

A=(~) (cosh (3) Nl (cosh a) Nv . (17) 
It is very easy to show that 

E CT «.=±1 E ff i, .=±1 ( X + ^ a i,j (T i',j){ 1 + VijVi+lj) 

(1 + y/zaf± ltj <Ti+ij) 

= 2 2 (l + zaija i+ltj ). (18) 

Hence Z and Z are equal. 

Now let us define the following set of constants {Oij}. 

0i,i = l 

n 

#l,n+l = E Q\,k + 1 
k=l 

m— 1 TV n 

E E + E + L ( 19 ) 
j=i fe=i k=i 



'm,n+l 



The solutions of these equations are 

e m)n = 2 { - m -^ N + n - 1 . (20) 

Now let us consider the following linear combination 

JV TV 

E E "<-Aj «m = o,±i. (21) 
»=i j=i 

We are going to show that (|21jl holds only if the integer coefficients 
rtjj are zero. It is very easy to show that 

TV N-1 

0N,N > I E E n id i,j\ V «M = °> ±1- ( 22 ) 
i=l i=l 

Hence if (|2~Tj) holds then n/v,jv must be zero. 
In the same way we can show that 

TV TV-2 

#TV,TV-1 > I E E n i,3°i,j\ Vn M = °' ±L ( 23 ) 
i=l j=l 
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Then 

tin.n— i must be zero. 
In general we can show that 

m—1 N n-1 

^2 n m ,j9ij\ Vriij = 0, ±1. (24) 
i=i j=i j=i 

Then n n>m must be zero. 

Hence if the linear combination defined in (|21[) is zero then all the 
coefficients rtj j must be zero. 

Now let us define the partition function: 

1 rTfj rTy _ 

Z' = A lim lim — — / dun / du>v Tl 
^T H -+ooT v -+ooT H T V Jo H Jo ' 11 

(l + y/zaija^J (l + x/zcJijaPj 

(1 + /iCTij) 

(l + -rr,'^ ) (l + e-^^^og) , (25) 

Let us consider periodic boundary conditions. In this case the bound- 
ary conditions of the constants {Oij} are given by 

9o,j = Onj 

Qi,0 = 6i,N- (26) 

We are going to show that this partition function is equal to Z. 
Let us remark that 

JJ(1 + Oi) = 1 + 5Z a ' +^2 a i a j + OjOj-afc H . (27) 

j i i<j i<j<k 

We are going to use this formula in order to expand the products: 

Ui,j (l + < '~'"" ■ (l + e^Wffj,] 

(l + e^^'o-^-) (l + e-^^aP^j . (28) 

Hence (|28|) is given by a sum of terms like 
. . . + e ^(E • , TJ TJ m '-i^.') n n>£ + ' ( 29 ) 
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where indices J and K can take the values 

J = R,L 

K = U,D. (30) 
Let us remark that n^j and m^j can take the values 

riij = 0, ±1 

tthj = 0,±1. (31) 

The primes over the sum and product symbols in (|29[) means that the 
indices do not take all their values. 

For instance, the partition function Z' depends on 9i j through the 
factor 

= (l + - + + a^ L +1J ) . (32) 

Let us remark that first and fourth term in the right hand side of the 
above equation correspond to riij = 0. The second term is defined by 
riij = 1 and the third term is related with Ujj = — 1 

Now we are going to perform the integration over ujh an d ujy and 
to take the limits in T# and Ty. The only terms of the expansion 
Q29|) that survive to the limits in Th and Ty are those with all the 
coefficients n^j and m^j equal to zero. Therefore factors like (|32|) 
becomes 

(l + ^ L +1J ). (33) 

Hence is also given by: 



'l + y/zOijO%J) (l + y/zatjoPj) (1 + fto^ 



1 + oSof+ij) (l + • (34) 

This is the partition function Z given in (jlfij) . 

3. Integral representation of the partition function 

We have shown that the partition function of the 2D Ising model 
with a magnetic field is equivalent to the partition function Z' given 
in l|25j>. Hence Z can be written as: 

I fTjj fTy 

Z = A\m\ lim duj H dwy ^ JJiQj, (35) 

T h -kx> Ty-voc ±h±V JO JO <rij 
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where Kij is given by: 

Kij = (l + s/za i:j afy (l + s[zoija\j 

(l + sfzoija^ [\ + y/za i:j afy (1 + ha i:j ) 

(l + e iuJHe ^o*) (l + e-^^a^) 

(i + e **i "' ■n'j) (l + e-^^-^g) . (36) 

Let us remark that Kij depends only on the spin variables defined at 
the vertex (i, j). Hence 

£n^=n£^ ( 37 ) 

We can perform the integration over the spin variables and the parti- 
tion function Z is given by 

1 rT H rT v _ 

Z = A]hn lim — — / du H du v [[Kij, (38) 
t h ^oo t v ^oo ±hJ-v Jo Jo 

where 

= £ £ £ £ £ 

(l + sTz°ij(jfj) (l + ^a id af :j ) 

(l + y/zvijO-Yj) (l + s/za i:j afy (1 + /wjij) 

(1 + e ^^.^5) (l + e-™* '- 1 *^ 

(l + e^^a^) (l + e-^^-^g) . (39) 

It is very easy to show that Kjj is given by 

Kij = 2 5 k { j, (40) 

where 

_|_ z "i e idi,j{^H+^v) _|_ z 2g— i(wfrSi_i j j+a;v^i,j-i) 
_|_ ^4 e i (u H 6»i j : +u v 6»i j- -cjjf 0j _ i j ; - w v 6i j • _ i ) 
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_|_ Z(j j IgtWySi,] _|_ e '^H"i,j _|_ g— lwv"i,j-l _|_ e —1^>B"i-l,j 

+ z 3 we' 
+ z 3 we* 
+ z 3 uje l 

+ z 3 we i ' WHS,J+ ' Jv ' 9i ' r ' J, ' 9,J '" 1 '. (41) 
Hence the partition function Z has the following representation: 

Z = 2 5N2 A lim lim -J— ^ du v exp 

Th^ooTv^oo Ih^V JO JO 

(42) 

where 

K^lni^ (43) 

and i^jj- is given by (|4Tj) 

Hence we have transform a sum over an arbitrary number of vari- 
ables into and integral over two variables, a double sum over two 
indices and two limits. 

4. Conclusions 

We have shown that the partition function of the 2D Ising model 
coupled to an external magnetic field can be represented by an inte- 
gral over a finite number of degrees of freedom. This result is very 
interesting because the 2D Ising model with magnetic field cannot be 
solved exactly. 

The remaining finite integration must be performed numerically. 
But it is difficult to perform numerical calculation with this repre- 
sentation of the partition function. The main problem is the total 
memory that we need in order to store the coupling constants Oij. If 
6ij are stored as 32 bits integers then 

- 2 31 < 9 itj < 2 31 . (44) 

If they are stored as 64 bits floating point then the maximum value of 



Vij is 

In 2D this bound means that 



< 2 1000 . (45) 



iV < 30 (46) 

but in 3D 

N < 10 (47) 
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